The Abundance of Atomic Sulfur in the Atmosphere of Io 



Lori M. Feaga 1 , Melissa A. McGrath 2 , and Paul D. Feldman 



ABSTRACT 

Observations with the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope 
have been used to constrain the atomic sulfur column density in Io's atmosphere. The S I A1479 dipole 
allowed and forbidden transition multiplets have been resolved for the first time at Io, enabling the study 
of both optically thick and thin transitions from a single atomic species. The allowed transitions con- 
tribute 62 ± 8% and the forbidden transitions 38 ± 8%, on average, to the total signal of the S I A 1479 
multiplets. Using the optically thick and thin transitions of S I A 1479 observed near the limbs of Io, we 
derive a tangential atmospheric sulfur column abundance of 3.6 x 10 12 cm -2 < M s < 1.7 X 10 13 cm -2 , 
which is independent of electron temperature and density. A low density SO2 atmosphere, M S02 ~ 5- 
lOx 10 15 cm~ 2 , consistent with that inferred from other recent observations, is most consistent with these 
bounds. 

Subject headings: planets and satellites: individual (Io) — ultraviolet: solar system — atomic data 



1. Introduction 

Io's atmosphere plays a key role in populating the 
jovian magnetosphere with plasma and replenishing 
the plasma torus with sulfur and oxygen (Spencer & 
Schneider 1996). This tenuous atmosphere, first de- 
tected by the Voyager 1 IRIS experiment over 20 years 
ago (Pearl et al. 1979), has proven to be a challenge 
to characterize quantitatively. It is composed primar- 
ily of S0 2 (Lellouch et al. 1990, 1992) and SO (Lel- 
louch et al. 1996). The atomic by-products of SO2, 
sulfur and oxygen, were first identified in the Io plasma 
torus (Brown 1981; Durrance, Feldman, & Weaver 
1983), and more recently detected in Io's emission line 
spectrum (Ballester et al. 1987). Sodium (Burger & 
Schneider 1998; Geissler et al. 1999; Bouchez, Brown, 
& Schneider 2000) and potassium (Brown 2001) are 
also present at much lower abundance than sulfur and 
oxygen. The observed emissions are produced when 
the atmospheric sulfur and oxygen are collisionally ex- 
cited by impinging electrons from the plasma torus. 
The emissions are therefore diagnostic of the plasma 
interaction with Io's atmosphere, providing quantita- 
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tive information about the electron temperature and 
density, in addition to the abundances of sulfur and 
oxygen. 

Despite the recent advances in quantifying Io's at- 
mosphere, including investigating the inhomogeneity 
and impact of volcanic activity (Feldman et al. 2000; 
McGrath et al. 2000; Spencer et al. 2000), analyz- 
ing the airglow (Geissler et al. 1999), and resolving 
the morphology (Roesler et al. 1999; Retherford et al. 
2000), and the fact that several sets of sulfur and oxy- 
gen emission line data have been acquired and ana- 
lyzed (Ballester 1989; McGrath et al. 2000; Oliversen 
et al. 2001; Wolven et al. 2001), the abundances are 
still not well constrained. Various observational mea- 
surements and theoretical models have placed loose 
bounds on the sulfur column density (M s ). Previ- 
ously, electron densities and temperatures were as- 
sumed in order to estimate M s , resulting in limits of 
2.2xl0 12 cm~ 2 < M s < 7xl0 15 cm" 2 (Ballester et 
al. 1987; McGrath et al. 2000). Recent observations 
make both optically thick and optically thin regimes 
available in a single exposure, allowing us to constrain 
M s without assuming a specific T e or n e . We present 
such a temperature and density independent method 
for determining the column density in this paper. 

Sulfur and oxygen abundances also provide a little- 
utilized diagnostic of the origin and nature of the SO2 
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atmosphere. Significant progress has been made in 
characterizing both the spatial distribution and tem- 
poral variability of the SO2 atmosphere (Ballester et 
al. 1994; Trafton et al. 1996; McGrath et al. 2000; 
Feldman et al. 2000; Spencer et al. 2000; Strobel 
& Wolven 2001). The current picture includes a 
higher column density equatorial, sub-solar region, 
Nso-2 ~ l-5xl0 16 cm~ 2 , with a strong poleward gra- 
dient, where Af S02 decreases to < 10 15 cm -2 (Feldman 
et al. 2000; Strobel & Wolven 2001) and column den- 
sities > 10 16 cm~ 2 above active volcanic regions (Mc- 
Grath et al. 2000; Spencer et al. 2000). It is obvious 
that both active volcanism and SO2 frost sublimation 
are important atmospheric sources, although the rela- 
tive importance of each mechanism has not been de- 
termined (Summers & Strobel 1996; Wong & Smyth 
2000; Moses, Zolotov, & Fegley 2001). The models 
show that sulfur can serve as a proxy for SO2, with 
enhanced abundance relative to oxygen during and 
shortly after volcanic activity as compared to a subli- 
mation atmosphere. Sulfur displays similar trends as, 
and scales well with, the SO2 abundance in most of 
the day side models. Determining and monitoring the 
sulfur abundance, can, in principle, provide a much- 
needed measure of the variability of Io's atmosphere. 

We describe the observations and explain the data 
reduction in Section 2. We then give the motivation 
for and present the analysis of the opacity of the sul- 
fur transitions and the atomic sulfur column density in 
Section 3. Finally, we discuss the implications of our 
findings in Section 4. 

2. Observations and Data Reduction 

The Space Telescope Imaging Spectrograph (STIS) 
aboard the Hubble Space Telescope (HST) has made 
it possible to study the sulfur emission near 1479 A in 
detail. The allowed multiplet, S I A1479 (3s 2 3p 3 4s' 3 D' 
3s 2 3p 4 3 P), and the intercombination multiplet, S I] 
A1479 (3s 2 3p 3 3d 5 D° -> 3s 2 3p 4 3 P), are both present 
in the data and are resolved for the first time. The 
analysis that follows focuses on the separation of the 
S I A 1479 multiple ts and its implications. Additional 
data collected with the International Ultraviolet Ex- 
plorer (IUE), from 1986-1988, and the Faint Object 
Spectrograph (FOS) on HST, from 1992-1996, are re- 
examined in light of the STIS data. We present only 
one representative spectrum per instrument in this pa- 
per. We concentrate specifically on the following sul- 
fur emission multiplets: the intercombination and al- 



lowed multiplets at S I A1479, the allowed multiplet 
S I A1814 (3s 2 3p 3 4s 3 S° -^3s 2 3p 43 P), and the in- 
tercombination doublet S I] A1900 (3s 2 3p 3 4s 5 S° -> 
3 s 2 3p 4 3 P). For all but the 1479 A multiplet, there are 
several years of data that confirm the trends and ratios 
found in our analysis. Table 1 summarizes the data 
analyzed, with the observations presented in this paper 
shown in bold. The morphology of the STIS obser- 
vations is identical to that described in Retherford et 
al. (2000), with resolved bright auroral emission re- 
gions near the Io equator at both sub- and anti-Jupiter 
locations. The equatorial spots are not resolved in the 
other observations presented here. 

The primary data set we analyze was acquired with 
STIS on HST using the 52" x 2" slit and the G140M 
grating, which produces a monochromatic image of 
Io at each emission line detected in the spectral range 
1440-1500 A. The set contains four ~19 minute expo- 
sures taken on 1998 August 21-22 while Io was near 
western elongation, Io orbital longitude 285°-305°. 
The coadded flux spectrum of the first two exposures is 
presented in Figure 1 . The data are acquired at a dis- 
persion of 0.05 1 A/pixel and spectral resolution of <~ 
0.5 A, allowing 6 blended components of the S I A 1479 
multiplets to be resolved. The detected features corre- 
spond to both the electric dipole allowed and forbid- 
den transitions shown in the Grotrian diagram in Fig- 
ure 2. Atomic data for the transitions are given in Ta- 
ble 2. The blended allowed transition wavelengths are 
1474.1 A, 1483.1 A, and 1487.2 A, while the blended 
forbidden transitions are at 1473.0 A, 148 1.7 A, and 
1485.6 A. 

These data were reduced in the following man- 
ner. The rectified (x2d) images were boxcar smoothed 
over a 3 x 3 pixel region and background subtracted. 
The wavelength dependent background level was de- 
termined by averaging the signal in a region spatially 
one Io diameter away from and comparable in size to 
the bright spectral features. Next, several rectangular 
sub-images, centered on the brightest emission areas 
of the anti- and sub-jovian equatorial features, were 
created in order to isolate these spots and to ensure no 
overlap of emission. Lastly, the data in the sub-images 
were compressed spatially into a one-dimensional ar- 
ray of flux vs. wavelength. Because the spots are not 
located at the center of the aperture, their flux spec- 
tra have systematic wavelength offsets, which we have 
corrected for. A Gaussian profile was fitted to each of 
the 6 resolved emission features to determine their rel- 
ative strengths. These Gaussians are overplotted on the 
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data in Figure 3. 

We re-analyze a set of IUE observations of Io which 
span the spectral range 1150-1950 A (Ballester 1989). 
Low resolution ^5 A exposures of 13-14 hours in 
length were made over a three year period with the 
short-wavelength (SWP) camera. We present a rep- 
resentative spectrum in Figure 4a which is centered 
on Io's western elongation on 1987 July 31; it was 
reduced following the same method as Ballester and 
smoothed using a 5 -element boxcar routine. An ad- 
ditional step was taken to remove the reflected light 
continuum longward of 1700 A by subtracting a solar 
spectrum scaled to the intensity of the underlying con- 
tinuum. 

The last set of Io observations we re-analyze was 
made with the FOS. These observations cover 1600- 
2300 A and were obtained over a period of four years 
using the G190H grating, which has a dispersion of 
1.45 A/diode. Exposures varied in length between 10 
and 20 minutes. The representative spectrum shown 
in Figure 4b was acquired while Io was emerging into 
sunlight from Jupiter's shadow on 1996 October 2. In 
order to minimize scattered light contamination, only 
the first 7 groups of the data, corresponding to Io in 
eclipse, were used. The data were reduced in the same 
manner as presented in McGrath et al. (2000) and were 
rebinned by 2 pixels and smoothed by 3. 

Several factors were considered in the analysis, 
including the influence of the emission morphology 
on the results, and possible contamination from SO2 
absorption. The emission morphology (e.g., one or 
two spots visible) depends on Io's orbital geometry 
(Roesler et al. 1999; Retherford et al. 2000), which 
may affect the absolute magnitude of the measured 
fluxes. When Io is at western elongation, as it is for 
part of the IUE and all of the STIS integration time, 
both equatorial spots are in view for at least part of 
the exposure. For the FOS egress data, when Io is in 
shadow emerging from behind Jupiter, only one equa- 
torial spot would be visible if resolved. We assume that 
the opacity and the S I A1479 analysis we present are 
independent of the emission morphology. Concerning 
contamination due to SO2 absorption, since the SO2 
is more concentrated at lower altitude than the atomic 
sulfur in the atmosphere, the absorption should be min- 
imal. Also, at the relevant wavelengths, the average 
SO2 transmission is close to unity. Therefore, the rel- 
ative line strengths of the sulfur multiplets should not 
be significantly affected. 



3. Analysis 

3.1. Allowed and Forbidden Components of S I 
A1479 

It is desirable to separate the allowed from the for- 
bidden transitions since the transition probabilities and 
the oscillator strengths for the two transitions are very 
different. The forbidden transitions are much more 
likely to be optically thin than the allowed ones, mean- 
ing two regimes of opacity can potentially be sampled. 
A quantitative determination of the relative contribu- 
tions of the allowed and forbidden transitions can also 
be used to help analyze numerous data sets for which 
the multiplets are unresolved. 

A synthetic sulfur emission spectrum is created for 
the STIS data set by adding together 13 Gaussians of 
the same FWHM, one at each of the central wave- 
lengths of the components of the S I A 1479 multiplet. 
Each Gaussian is weighted by its optically thin relative 
strength as reported in Morton (1991) and then the en- 
tire synthetic spectrum is scaled to the intensity level 
of the data. A fit of the synthetic spectrum to the data 
is shown in Figure 3. The flux from each line is de- 
termined by integrating the area under the Gaussian 
curves in Figure 3 and is reported in Table 3. The al- 
lowed transitions contribute an average of 62 ± 8% 
and the forbidden transitions an average of 38 ± 8% to 
the total signal of the S I A1479 emission. Although 
the anti- and sub-jovian equatorial spots differ in flux 
by up to a factor of 2, the relative line ratios and the 
contribution of the allowed and forbidden components 
are consistent throughout the data set. 

3.2. Opacity of the Observed S I Emissions 

A synthetic sulfur emission spectrum is also created 
for the IUE and FOS data sets. A unique FWHM is 
chosen for each data set to match the resolution of the 
instrument. Again, Gaussians are placed at the central 
wavelengths for all lines of a given sulfur multiplet, 
weighted by its optically thin relative strength (Morton 
1991), and scaled to the data. If the data match the 
synthetic spectrum, sulfur is considered optically thin 
in that emission; otherwise, it is considered optically 
thick. 

The fits for S I A1814 and S I] A1900 are over- 
plotted on the data in Figure 4. If optically thin, the 
transitions of the S I A1814 multiplet should be in a 
5.2 : 3 : 1 ratio and the S I] A 1900 doublet will be in a 
3.6 : 1 ratio. The ratios are calculated from their transi- 
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tion probabilities listed in Table 2. We observe a ratio 
of 3.4: 5.6:1 and 4.5:5:1 for the 1814 multiplet us- 
ing the IUE and FOS data respectively, and 3 : 1 and 
3.4 : 1 in the 1900 doublet. Our results therefore indi- 
cate that the S I] A 1900 doublet is optically thin and 
the S I A1814 multiplet is optically thick. Moreover, 
we have analyzed all the IUE and FOS data listed in 
Table 1 and find that in all existing observations of Io, 
the 1900 A doublet is observed to be optically thin and 
the 1814 A multiplet is observed to be optically thick. 

The allowed transitions of the S I A 1479 multiplet 
have a theoretical ratio of 5.2: 3.1 : 1 when optically 
thin (Morton 1991). We observe a ratio of 5.3 : 3.8 : 1 
in the STIS data for the S I A1479 allowed multi- 
plet. Initially, our line ratios seemed to indicate that 
the S I A 1479 allowed multiplet is optically thin since 
the blended line ratios are consistent with the optically 
thin values; however, this contradicts the fact that the 
S I A1814 multiplet, which has the same ground state 
and similar oscillator strength as S I A1479, is con- 
sistently optically thick in all existing observations of 
Io. A more detailed analysis of the S I A1479 allowed 
transitions shows that because the lines are blends at 
STIS resolution, optical depth effects can be hidden 
by the blending. The shorter wavelength transitions of 
a given blend, for example 1473.99 A, 1474.38 A, and 
1483.04 A, are preferentially depleted at high optical 
depth. The flux in these transitions is redistributed in 
the multiplet such that the fluxes in the longer wave- 
length transitions of the blend, 1474.57 A and 1483.23 
A, are enhanced, and the 5 : 3 : 1 ratio of the blends is 
unaffected by this flux redistribution. 

3.3. Constraints on M s 

The optical depth of a transition is directly propor- 
tional to both the column density and the absorption 
cross section, r = TVer, with r = 1 serving as the 
conventional boundary between the optically thin and 
thick regimes. Therefore, generally speaking, an upper 
limit can be placed on the sulfur column density in the 
optically thin regime, and a lower limit in the optically 
thick regime. If a is known, the sulfur column density 
can be constrained without assuming a specific T e or 
n e . 

The absorption cross section of photons in the core 
of an optically thin transition where only thermal 
Doppler broadening and radiative damping are consid- 
ered is directly proportional to the oscillator strength, 
fij, of the transition (Osterbrock 1989). Assuming t 



= 1, 

r 1 m e cAi i D rrieC [2kT 

N =-= = 4 / ^v . (1) 

a V ne Jij V ne Jij V m s c z 

Several consistent values of the oscillator strength are 
published for the dipole allowed transitions S I A1814 
and S I A1479. Only one value of the oscillator 
strength exists for the majority of the forbidden tran- 
sitions of S I] A 1479. For the S I] A 1900 transition, 
only two values of the oscillator strengths and transi- 
tion probabilities exist (Muller 1968; Tayal 1998-see 
the discussion in McGrath et al. (2000) on the order of 
magnitude disagreement between these values). The 
values and references for the transition probabilities 
and oscillator strengths we use are given in Table 2. 

Using equation (1) with T = 1000K and the oscil- 
lator strength of fa = 0.090 (Morton 1991) for the S I 
A 1479 allowed multiplet, the lower bound for the sul- 
fur column density is3.6xl0 12 cm~ 2 . A similar lower 
bound of 2.8 x 10 12 cm~ 2 is found using the value 
of 0.093 (Tayal 1998) for the S I A1814 multiplet. An 
upper bound of 1.7x 10 13 cm~ 2 is found using the os- 
cillator strength of the dominant S I] A 1479 forbid- 
den transition line, 1472.97 A, listed in Morton (1991) 
as fij = 0.019. This result is lower than the previous 
upper bound for the atomic sulfur column density at 
Io, 7x 10 15 cm -2 , by more than 2 orders of magnitude 
(McGrath et al. 2000). An explanation of this discrep- 
ancy is given in the discussion. 

3.4. Constraints on T e and n e 

Knowledge about the sulfur column density in turn 
allows us to place constraints on the electron temper- 
ature and density of the impacting plasma. The ob- 
served line of sight brightness is given by 

B i:j = J n s n e Qy(T e )d/ (2) 

where n s and n e are the sulfur and electron densities, 
respectively, and Qij(T e ) is the temperature depen- 
dent electron excitation rate coefficient. Although it 
has been shown by Combi et al. (1998) and Saur et al. 
(1999) that the plasma flow and interaction around Io 
creates varying electron density and temperature pro- 
files at Io, we assume that n e and T e are constant along 
the line of sight, invert equation (2), and arrive at an ap- 
proximate relationship between the sulfur column den- 
sity and the measured emission line brightness, 
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where Bij is in units of Rayleighs, n e is in units of 
cm -3 , and Qij is in units of cm 3 s -1 . For a particular 
transition from an upper state j to a lower state i, the 
rate coefficient is (Osterbrock 1989) 



8.63 x 10 -6 Oye- £ «/ feTe 

(°Cij\-L e) 



(4) 



where fiy is the thermally averaged collision strength, 
uji is the statistical weight of the lower state, is 
the transition energy in units of eV, k is the Boltzmann 
constant, and T e is the electron temperature in units of 
Kelvin. These equations apply only to optically thin 
transitions, therefore, the S I] A1479 forbidden multi- 
plet or the S I] A 1900 intercombination doublet must 
be used in determining the electron temperature and 
density. As there are no electron impact excitation 
cross-sections or collision strengths available for the 
S I] A1479 forbidden multiplet, the S I] A1900 inter- 
combination doublet is used for the analysis. We use 
the electron excitation rate presented in Figure 9a of 
McGrath et al. (2000). 

Combining the sulfur column density bounds from 
the S I A 1479 analysis, the line brightnesses of the 
S I] A 1900 intercombination doublet from the IUE 
data when Io was west of Jupiter, and equation (3), 
the density and temperature of the impacting electrons 
can be estimated. For an electron temperature range of 
0.8-100 eV and several impacting electron density val- 
ues between 500 and 10000 cm -3 , the column density 
is plotted in Figure 5. Both a low temperature solution, 
1.5 eV < T e < 7 eV, and a high temperature solution, 
7 eV < T e < 100 eV, are possible. The canonical val- 
ues of 5 eV and 2000 cm -3 are consistent with the low 
temperature solution. Assuming the low temperature 
solution implies 1000 cm -3 < n e < 5000 cm -3 . 

4. Discussion 

The STIS observations presented here and used to 
calculate the atomic sulfur column density are tan- 
gential cuts through Io's atmosphere centered on the 
equatorial spots; therefore, the column densities are 
also tangential. Wolven et al. (2001) has recently pub- 
lished a collection of four years of HST/STIS low res- 
olution observations of near Io emissions and extended 
emissions in the 1150-1730 A wavelength range. The 
photon fluxes measured at 1 Io radius for the S I 
A1479 multiplets varied between 3 x 10 -4 cm -2 sec -1 
and 3x 10 -3 cm -2 sec -1 , corresponding to eclipse 
and western elongation, respectively. These values 



are consistent with our HST/STIS photon fluxes for 
S I A1479, implying that the tangential sulfur col- 
umn that Wolven et al. (2001) were sampling should 
have values for M s similar to ours. McGrath et al. 
(2000) determined the vertical column density of S, 
SO, and SO2 above specific locations on Io, including 
the Pele volcano. With an estimated collision strength 
and electron impact excitation rate coefficient for the 
S I] A 1900 doublet, and assuming canonical plasma 
torus values of electron density and temperature, 1000 
cm -3 and 5 eV, they determine an M s of 1 x 10 14 cm -2 
at Pele, which is much larger than our upper bound of 
1.3xl0 13 cm -2 . A column density that is significantly 
larger than our value should result in additional opti- 
cal depth effects. In Figure 6, the Pele flux spectrum 
published in Figure 2 of McGrath et al. (2000) is com- 
pared with the optically thin FOS synthetic spectrum 
described earlier, and confirms that the S I A1814 mul- 
tiplet is more optically thick over Pele than in our IUE 
and FOS spectra, lending credence to the significantly 
higher value of M s over Pele despite the uncertainty 
in the S I] A 1900 electron excitation rate. This im- 
plies that M s is larger in volcanic plumes than in the 
equatorial spots. 

In the 1-D, steady state photochemical atmospheric 
models presented by Summers & Strobel (1996), both 
high and low density SO2 atmosphere sub-solar pro- 
files with various levels of vertical eddy mixing are 
calculated. We have extracted the vertical column den- 
sities of several atmospheric components from their 
published profiles. The A^ s values corresponding to the 
high density SO2 cases, M S02 = 10 18 cm -2 , were much 
higher than the limits from our analysis. For exam- 
ple, with a low vertical eddy mixing coefficient of k = 
10 6 cm 2 sec -1 , M s ~ 2x 10 17 cm -2 ; for a high vertical 
eddy mixing coefficient of k = 10 9 cm 2 sec -1 , ~ 
9x 10 14 cm -2 . On the other hand, for their low den- 
sity case (M S02 = 8 x 10 15 cm -2 with high eddy mixing 
k = 10 9 cm 2 sec -1 ) M s ~ 6x 10 12 cm -2 . The corre- 
sponding tangential column density, assuming spheri- 
cal symmetry and using the atmospheric sulfur num- 
ber density profile in Summers & Strobel (1996), is on 
the order of 4x 10 13 cm -2 . This value is much more 
consistent with our results and the general overall pic- 
ture of the SO2 atmosphere described in the introduc- 
tion. In a 3-D high density SO2 model atmosphere 
presented by Wong & Smyth (2000), they assume an 
initial sub-solar M S02 of 6x 10 17 cm -2 , which results 
in a tangential sulfur column density at the terminator, 
M s ~ 2x 10 15 cm -2 , inconsistent with our upper limit 
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to Af s . In the model in which the input SO2 column is 
reduced to a lower density value of 7x 10 16 cm~ 2 , the 
tangential J\f s is more consistent with our results at a 
value of6x 10 13 cm~ 2 . 

In order to predict what the atmospheric driving 
mechanism is, Moses, Zolotov, & Fegley (2001) have 
modeled both an SO2 frost sublimation atmosphere 
and a Pele-type volcanically driven atmosphere. Of 
the two models, our J\f s value is more consistent with 
the SO2 frost sublimation atmosphere simply because 
this model has a lower sulfur column density. How- 
ever, the true test to the Moses et al. (2001) models 
is to compare the sulfur to oxygen column abundances 
extracted from a single exposure. When the sulfur is 
more abundant than the oxygen, the models imply a 
volcanically driven atmosphere. None of the medium 
resolution data presented here had a simultaneous oxy- 
gen detection, so the ratio is not currently available, 
although this sulfur to oxygen relationship will be an- 
alyzed in future work with a more extensive data set. 
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Table 1 
Summary of Io Observations 



Observation 


UT Date 


UT Time 


Io Diameter 


Exposure Time 






(hh:mm(:ss)) 


(") 


(s) 


IUE 


swp28708 


1986 Jul 18 


21:49 


1.14 


48600 


swp2943() 


1986 Oct 12 


15:29 


1.21 


49800 


swp31440 


1987 Jul 31 


18:53 


1.06 


50400 


swp31447 


1987 Aug 1 


16:57 


1.07 


48000 


swp34341 


1988 Sep 27 


17:02 


1.12 


48600 


swp34343 


1988 Sep 28 


14:49 


1.12 


47700 


HST/FOS 


y0w00203-4 


1992 Mar 22 


13:18:10 


1.11 


1680 


y0w05402 


1992 May 16 


10:17:41 


0.97 


1245 


ylal0801 


1993 Aug 4 


02:44:55 


0.84 


851 


y2wta201 


1995 Oct 9 


12:48:19 


0.86 


738 


y3cwbl01-5 


1996 Aug 1 


07:20:06 


1.17 


4920 


y3eta201 


1996 Oct 2 


19:24:22 


0.98 


800 


HST/STIS 


o4xm03050 


1998 Aug 21 


22:07:17 


1.23 


1155 


o4xm03060 


1998 Aug 21 


22:30:02 


1.23 


1155 


o4xm03()7() 


1998 Aug 21 


23:42:42 


1.23 


1145 


o4xm03()8() 


1998 Aug 22 


00:05:17 


1.23 


1145 
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Table 2 
Atomic Data for SI 



Multiplet 






Actual A 


Transition Probability 


Oscillator Strength 


Reference 










(A) 


(s- 1 ) 






B D _ 


+ 3 P 


















3 


2 


1472.97 


4.20 X 10 7 


1.91 X 10~ 2 


1 






2 


2 


1473.02 


6.15 x 10 5 


2.00 x 10" 4 


2 






1 


2 


1473.07 


3.91 x 10 4 


7.63 x 10 -6 


2 






2 


1 


1481.66 


1.70 x 10 7 


9.33 x 10- 3 


1 






1 


1 


1481.71 


4.28 x 10 5 


1.41 x 10- 4 


2 









1 


1481.74 


3.84 X 10 3 


4.21 x 10~ 7 


2 






1 





1485.62 


2.30 X 10 6 


2.28 x 10 -3 


1 


3 D U _ 








1478.50 


1.65 x 10 8 


9.01 x 10" 2 


1 






3 


2 


1473.99 


1.60 x 10 8 


7.30 x 10- 2 


1 






2 


2 


1474.38 


5.00 x 10 7 


1.63 x 10 -2 


1 






1 


2 


1474.57 


6.20 x 10 6 


1.21 x 10- 3 


1 






2 


1 


1483.04 


1.20 x 10 8 


6.60 x 10~ 2 


1 






1 


1 


1483.23 


7.50 x 10 7 


2.47 x 10~ 2 


1 






1 





1487.15 


8.70 X 10 7 


8.65 x 10 -2 


1 


3 S U _ 


+ 3 P 






1813.73 


5.67 x 10 8 


9.33 x 10" 2 


2 






1 


2 


1807.31 


3.17 X 10 8 


9.32 x 10 -2 


2 






1 


1 


1820.34 


1.88 X 10 8 


9.33 x 10- 2 


2 






1 





1826.24 


6.23 X 10 7 


9.34 x 10" 2 


2 


b S u _ 




















2 


2 


1900.29 


6.60 X 10 4 


3.60 x 10 _B 


3 






2 


1 


1914.70 


1.80 X 10 4 


1.70 x 10- 5 


3 



REFERENCES.— (1) Morton 1991; (2) Tayal 1998; (3) Muller 1968. 
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Table 3 
Measured Fluxes for SI 



Multiplet 


Central A 

(A) 


STIS Flux* 

(io- 3 


IUE Flux* 
photons cm~ 2 


FOS Flux* 


BpO _> 3p 


1479 

1473.0 
1481.7 
1485.6 


0.51 ± 0.09 

0.376 ± 0.083 
0.107 ±0.027 
0.026 ± 0.007 






a D u — ► 3 P 


1479 

1474.1 
1483.1 
1487.2 


0.76 ± 0.11 

0.394 ± 0.059 
0.289 ± 0.087 
0.075 ± 0.019 






3 S U _> 3p 


1814 

1807.3 
1820.3 
1826.2 




5.4 ± 0.9 

1.7 ±0.4 

2.8 ± 0.4 
0.5 ± 0.2 


1.1 ± 0.2 

0.45 ± 0.08 
0.50 ± 0.07 
0.10 ±0.03 


bgl) _,. 3p 


1900 

1900.3 
1914.7 




5.1 ± 0.9 

3.8 ±0.5 
1.3 ±0.4 


2.0 ± 0.2 

1.48 ±0.09 
0.43 ± 0.08 



*Flux values are from entries in bold in Table 1. 
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1480 1485 
Wovelength (A) 



Fig. 1 . — Data extracted from the sum of the first two 
HST/STIS G140M exposures in the near vicinity of Io. 
The multiplet components of the S I A 1479 emission 
are seen. The brightest emission is located at the anti- 
and sub-jovian equatorial spots (top) with the compass 
showing the direction of Jupiter and Io North. The flux 
spectrum corresponding to the sub-jovian spot is plot- 
ted (bottom) with the data in black and the propagated 
errors in red. 



7 -179 ^ F 1-1:11 I E 



51 1479 A ALLOWED 



-6- 

7|7 



'•0 2 
3 Dl 



'F'O 

7 I ' 



Fig. 2. — Grotrian diagram for the S I A1479 dipole 
allowed and forbidden multiplets. All transition lines 
are shown. The lines indicated as blends are not indi- 
vidually resolved in the STIS data. 




1480 1485 
Wavelength (A) 



Fig. 3. — Synthetic fit to the flux spectrum of the STIS 
data previously shown in Fig. 1. The forbidden mul- 
tiplet fit is in red, the allowed in blue, and the total 
fit in boldface black. The resolved transition lines and 
blends are indicated. 




1800 1850 1900 

Wavelength (A) 



950 



Fig. 4. — Flux spectra of Io's atmospheric emissions in 
the wavelength range 1750-1950 A taken with (a) IUE 
and (b) HST/FOS. S I A1814 and S I] A 1900 emissions 
are present and are overplotted with an optically thin 
profile of the multiplets in bold. The optically thin pro- 
file is normalized to the blend of the 1820 A and the 
1826 A line. This comparison establishes that the S I] 
A1900 transition is optically thin whereas S I A1814 is 
optically thick. 
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Fig. 5. — Several sulfur column density curves cor- 
responding to different values of the electron density 
are plotted for the average measured brightness of S I] 
A 1900 taken from the IUE data with Io west of Jupiter. 
The upper and lower limits of the sulfur column den- 
sity from our analysis of the STIS S I A 1479 data, 
3.6xl0 12 cm~ 2 < M s < 1.7xl0 13 cm- 2 , are in- 
dicated by the dashed lines. The canonical T e and 
n e values of 5 eV and 2000 cm~ 3 give an Af s of 
5x 10 12 cm~ 2 , consistent with our bounds. 



1996 August 1 HST/FOS Data (Pele) 




1750 1800 1850 1900 1950 

Wavelength (A) 



Fig. 6. — Pele volcano flux spectrum from McGrath et 
al. (2000). A synthetic fit for optically thin neutral sul- 
fur emission, used to analyze the FOS data presented 
earlier in this paper, is overplotted on the data and is 
normalized to the 1826 A line. This illustrates that the 
S I A1814 multiplet experiences severe optical thick- 
ness effects while the S I] A 1900 doublet does not. 
This indicates that the sulfur column density is larger 
over Pele than the equatorial spots. 
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